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ABSTRACT 


Wave  loads  predicted  by  the  extended  SHIPMO  computer  program  were 
compared  with  experimental  data  in  an  attempt  to  assess  the  range  of 
validity  and  to  estimate  the  uncertainty  of  the  predictions.  In  general, 
the  agreement  between  theory  and  experiment  was  good  for  the  vertical 
bending  moment,  but  fair  to  poor  for  the  horizontal  bending  moment.  The 
predicted  torsional  moment  was  found  to  be  in  gross  disagreement  with  the 
measured  values  in  terms  of  both  magnitude  and  trend.  It  was  not  possible 
to  find  any  correlation  of  the  pattern  of  the  discrepancy  between  theory 
and  experiment  with  the  principal  design  characteristics,  such  as  hull 
form  or  weight  distribution.  The  present  comparison  shows  that  the 
prediction  of  ship  motions  in  oblique  seas  needs  to  be  further  improved. 
Also,  systematic  model  tests  covering  a  wide  range  of  the  variables  are 
required  to  assess  the  uncertainty  of  the  measurement,  and  hence,  to 
determine  the  range  of  confidence  of  the  predicted  loads  more  clearly. 


RESUME 


On  a  compare  les  charges  exercOes  par  les  vagues,  telles  que 
prOvues  par  le  programme  informatique  augmente  SHIPMO,  aux  donndes 
exp&rimentales ,  dans  un  effort  pour  Ovaluer  la  gamme  de  validity  et 
estimer  1' incertitude  des  provisions.  En  gOnOral,  1 'accord  entre  la 
theorie  et  l'experience  est  bon  pour  le  moment  flOchissant  vertical,  mais 
de  passable  A  mOdiocre  pour  le  moment  horizontal.  Le  moment  de  torsion 
pr&vu  s'&carte  beaucoup  des  mesures,  tant  en  amplitude  qu'en  direction. 

On  n'a  pu  Otablir  de  correlation,  quant  it  la  tendance  des  Ocarts  entre  la 
theorie  et  l'experience,  en  ce  qui  a  trait  aux  caractOristiques 
principales  de  conception  comme  la  forme  de  la  coque  ou  la  repartition  du 
poids.  La  pr&sente  comparaison  revfele  qu'une  amelioration  des  provisions 
des  mouvements  d'un  navire  qui  prend  les  vagues  de  c6t0  s' impose.  11  faut 
aussi  proceder  A  des  essais  systematiques  sur  maquettes  qui  portent  sur 
une  vaste  gamme  des  variables  pour  evaluer  1' incertitude  des  mesures  et 
determiner  ainsi  avec  plus  de  precision  l'intervalle  de  confiance  des 
charges  prOvues. 
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NOTATION 


aft  perpendicular 


R. 

Re{*} 


block  coefficient,  V/(LBT) 

bending  moment  coefficient  (■  VBM/pgBL2a,  or  HBM/pgBL2a) 
shearing  force  coefficient  (*  shearing  force/pgBIa) 
torsional  moment  coefficient  (“torsional  moment/pgBL2a) 
hydrodynamic  force  and  moment  on  portion  of  hull  due  to  body 
motion 

exciting  force  and  moment  on  portion  of  hull 
external  force  and  moment  on  portion  of  hull 
Froude  number 
forward  perpendicular 
horizontal  bending  moment 

inertia  force  and  moment  on  portion  of  hull 

length  between  perpendiculars 

restoring  force  and  moment  on  portion  of  hull 

real  part  of  {'} 

draft 

speed 

wave  load  component 
complex  amplitude  of  Vj 
real  part  of  VjQ 
imaginary  part  of  VjQ 
vertical  bending  moment 

complex  variable  for  either  of  Vj,  Ij,  Rj,  Ej,  Dj 

complex  amplitude  of  Xj 

real  part  of  XjQ 

imaginary  part  of  XjQ 

submerged  sectional  area 

two-dimensional  sectional  added  mass  coefficient 
sectional  beam  at  waterplane 
two-dimentional  sectional  damping  coefficient 
sectional  Froude-Kriloff  force 


acceleration  of  gravity 
sectional  diffraction  force 

/=r 

sectional  mass  moment  of  inertia  about  x-axis 

subscripts  indicating  modes  (j,k  =  1,2,. ..,6  refer  to  surge, 

sway,  heave,  roll,  pitch  and  yaw,  respectively) 

wave  number  ■  to2/g 
w 

sectional  mass  per  unit  length 
sectional  metacentric  height 
time 

coordinate  system  as  defined  in  Figure  1 
z-coordinate  of  sectional  centre  of  gravity 
wave  amplitude 

principal  sea  direction  relative  to  ship's  heading,  Figure  3 

displacement  volume 

motion  displacement 

(“  dn./dt)  velocity 
3 

(■  d2n j/dt2 )  acceleration 
phase  angle  of  wave  load 
wavelength 

variable  of  integration  in  x-direction 
mass  density  of  water 
frequency  of  wave  encounter 
absolute  wave  frequency. 
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INTRODUCTION 


One  of  DREA's  ship  dynamics  research  objectives  is  to  provide  the 
warship  designer  with  analytical  "tools"  which  enable  him  to  estimate  the 
motions  and  the  associated  wave  loads  on  ships  at  sea.  Recently, 
additional  subroutines1  to  calculate  the  wave  loads  have  been  developed 
and  integrated  into  the  ship  motion  computer  program  SHIPMO2.  The  new 
subroutines  calculate  the  root-mean-square  values  of  wave  loads  in  five 
degrees  of  freedom  in  irregular  short-crested  seas  as  well  as  the  transfer 
function  for  wave  loads  in  regular  seas. 

The  purpose  of  this  Technical  Memorandum  is  to  report  the  results 
of  a  study  which  compared  the  wave  loads  predicted  by  the  extended  version 
of  SHIPMO1  with  those  measured  in  experiments. 

The  method  of  analyses  of  motions  and  wave  loads  adopted  in 
SHIPMO  is  essentially  based  on  the  linear  strip  theory  due  to  Salvesen  et 
als.  However,  the  coefficients  in  the  equations  of  lateral 
anti-symmetric  motions,  including  roll,  are  modified  by  taking  into 
account  the  contribution  from  the  lifting  surfaces  on  the  hull  (stabilizer 
fins,  bilge  keels,  rudders,  propeller-shaft  brackets,  etc.),  the  details 
of  which  are  found  in  Reference  4.  The  summary  of  the  method  of  analysis 
of  wave  loads  in  regular  waves  is  presented  in  Section  2. 

In  Section  3,  the  predicted  wave  loads  are  compared  with  five 
sets  of  published  experimental  data. 


2.  METHOD  OF  ANALYSIS 


Figure  1  shows  the  right-handed  moving  cartesian  coordinate 
system  used  for  the  calculation  of  wave  loads  and  motions.  The  origin  is 
located  at  the  equilibrium  position  of  the  centre  of  gravity  of  the  ship. 
The  x-axis  points  to  the  bow  and  the  y-axis  extends  to  port,  both  being 
parallel  to  the  undisturbed  free  surface;  z  is  measured  vertically 
upward.  The  coordinate  system  maintains  a  fixed  orientation  with  respect 
to  the  undisturbed  free  surface  while  it  translates  with  the  mean  velocity 
of  the  ship.  The  dynamic  wave  load  components  acting  on  a  ship  cross 
section  are  defined  in  Figure  2. 

The  sea  direction  S6  is  defined  in  Figure  3  as  the  angle 
between  the  mean  velocity  vector  of  the  ship  and  the  principal  direction 
of  propagation  of  the  wave  train.  Thus,  $s  *  0°  for  following  seas, 

Bs  *  90°  for  beam  seas  from  port,  and  Bs  *  180°  for  head  seas. 
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Wave  Loads  in  Regular  Waves 


Consider  a  ship  moving  in  a  train  of  long-crested,  or  unidirectional 
regular  waves  of  a  particular  absolute  frequency,  ww,  on  a  straight 
course  at  a  fixed  orientation,  8S«  to  the  direction  of  the  waves” 
advance. 


Due  to  its  forward  motion,  the  ship  encounters  the  waves  with  a 
frequency  of  encounter  0)  which  is  related  to  the  absolute  wave  frequency 
via: 


where 


to  =  U)w  -  kw  UcosSs 

kw  =  wave  number  =  oo2/g 

w 

U  =  ship  speed 


(1) 


8S  *  sea  direction 

The  resulting  wave  loads  as  well  as  motions  will  then  be  harmonic 
functions  of  time  with  the  frequency  given  by  Equation  (1),  and  having  a 
solution  of  the  form: 

X j ( t)  =  Xj0exp(itot)  *  (XjR  +  iXjj)exp(iu)t)  (2) 

where  Xj  *  complex  variable  indicating  formal  solution  for  the 
motion  or  load  in  the  jth  mode 

XjQ  *  complex  amplitude  of  Xj 

Xjr»  XjI  =  real  and  imaginary  parts  of  XjQ,  respectively 
i  -  /=T 


to  *  frequency  of  wave  encounter 

t  *  t ime 


In  order  to  facilitate  the  algebraic  work,  the  computation  is 
performed  in  the  complex  field. 

In  expressions  involving  complex  numbers,  only  the  real  parts  are 
physically  realizable;  thus, 

Re{Xj  (t)}  =  XjRcoso)t  -  Xjjsinwt 

■  IXj0lcos((jjt  +  @j)  (3) 
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where  Re{Xj  (t)}  =  real  part  of  Xj 

2  2 

I X j 0 1  =  (Xjr  +  Xjj  ) l'2  =  amplitude 
6j  =  arctan  (Xjj/XjR)  =  phase  angle 

The  total  force  acting  on  a  cross  section  of  a  ship's  hull  may  be 
divided  into  two  components:  inertial  force  and  external  force.  Since 
inertial  forces,  the  magnitude  of  which  are  equal  to  the  mass  times 
acceleration,  act  in  the  directions  opposite  to  those  of  the  external 
forces,  the  net  loading  per  unit  length  of  a  ship's  hull  is  equal  to  the 
difference  between  these  two  forces: 

dVj/d£  =  dl j/d£  -  dFj/d£  (4) 

where  dVj/d£  =  net  load  in  jth  mode  per  unit  length, 

Ij/d£  =  local  inertial  force  in  jth  mode  per  unit  length, 

dFj/d£  ■  local  external  force  in  jth  mode  per  unit  length. 

The  modes  j  of  the  loading  are  defined  as  follows: 

j  “  1  axial  force 

j  *  2  horizontal  shearing  force 

j  =  3  vertical  shearing  force 

j  =  4  torsional  moment 

j  *  5  vertical  bending  moment 

j  =  6  horizontal  bending  moment 

The  net  wave  loads*,  Vj,  acting  on  an  arbitrary  cross  section 
are  obtained  by  integrating  the  load  per  unit  length  of  Equation  (4)  from 
one  end  of  the  hull  to  the  specified  cross  section. 

Thus  Vj  ■  / (dV j/d£)d£  =  net  wave  load 

Ij  =  /(dlj/d£)d£  =  inertial  forces  and  moments 

Fj  =  /(dFj/d£)d{j  =  external  forces  and  moments 


Hence,  from  equation  (1) 


!j  -Fj 


(j  =  1, . . . ,6) 


*  When  the  motion  of  an  entire  ship  is  considered,  the  external  forces  are 
exactly  balanced  by  the  inertia  forces,  so  that  the  sum  of  all  forces 
acting  on  the  ship  is  zero.  However,  the  two  kinds  of  forces  acting  on 
an  elemental  section  of  the  hull  are  generally  not  balanced. 
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It  can  therefore  be  seen  that  in  order  to  determine  the  wave 
loads  at  a  particular  cross  section,  the  forces  and  torsional  momemts 
acting  over  the  portion  of  the  hull  lying  on  one  side  of  the  cross  section 
in  question  must  be  known  &  priori,  and  that  unlike  the  calculation  of 
motions,  in  which  a  set  of  equations  must  be  solved  simultaneously ,  the 
calculation  of  wave  loads  on  a  cross  section  consists  simply  of  the 
summation  of  the  predetermined  forces  and  torsional  moments  acting  on  one 
side  of  the  section.  The  choice  of  the  side  in  any  case  is  of  course 
arbitrary.  In  the  subprograms,  the  summation  is  performed  over  the  hull 
forward  of  the  cross  section  under  consideration. 

The  axial  force  (j  =  1),  the  force  acting  in  the  direction  of  the 
x-axis,  is  assumed  negligible3  and  is  not  considered  herein. 

The  method  of  analysis  adopted  in  the  subprograms  is  the 
stripwise  computational  technique  of  Salvesen  et  al3.  The  forces  and 
moments  acting  on  cross  sections  are  derived  in  terms  of  motions 
(displacement,  velocities,  and  accelerations)  and  hydrodynamic 
coefficients  such  as  added-mass,  added  moment  of  inertia,  and  damping 
coefficients.  The  data  that  are  extraneous  to  the  calculation  of  motions, 
but  required  for  the  wave  load  calculation,  are  section  masses,  sectional 
radii  of  gyration  for  roll,  and  positions  of  sectional  centres  of  gravity. 

Outlined  below  is  the  computational  technique  implemented  in  the 
subprograms. 

The  hull  is  divided  into  a  number  of  equally  spaced  ordinate 
stations,  with  No.  0  station  at  the  forward  perpendicular  (FP).  The  total 
number  of  stations  may  be  up  to  21.  Each  hull  section,  or  "strip",  spans 
two  consecutive  half-stations  and  is  considered  to  be  a  cylinder  with  a 
constant  cross  section  identical  to  that  of  the  integral  station  at  its 
longitudinal  midpoint.  The  extreme  fore  and  aft  strips,  which  span  the 
integral  stations  at  the  perpendiculars  and  the  adjacent  half-stations, 
have  cross  sections  identical  to  those  at  the  perpendiculars. 

For  the  purpose  of  wave  load  calculation,  the  following 
assumptions  are  made: 

(i)  that  the  ship  is  a  rigid  body; 

(ii)  that  the  external  forces  and  moments  acting  on  each  strip  as 
well  as  the  mass  of  the  strip  are  concentrated  in  the  plane  of 
the  integral  station  contained  in  the  strip. 


The  external  forces  and  moments  are  caused  by  the  water  exerting 
forces  on  the  hull,  and  consist  of  three  components:  hydrostatic  restoring 
forces  and  moments  due  to  the  displacement  of  the  rigid  hull  from  the 
equilibrium  position,  exciting  forces  and  moments  due  to  the  oncoming 
waves,  and  the  hydrodynamic  forces  and  moments  due  to  the  ship's  motion. 
Hence , 

Fj  =  Rj  +  Ej  +  Dj  (j  =  2, . 6)  (6) 

where  Rj  =  hydrostatic  restoring  forces  and  moments 

Ej  =  exciting  forces  and  moments 

Dj  =  hydrodynamic  forces  and  moments  due  to  ship  motions 

From  Equations  (5)  and  (6),  the  wave  load  in  the  jth  mode  on  a 
cross  section  can  be  expressed  as: 

Vj  =  Ij  -  Rj  -  Ej  -  Dj  (j  =  2, . 6)  (7) 

The  following  expressions  for  the  I j ,  R j ,  E;  and  Dj  are 
given  in  Reference  3,  to  which  the  reader  is  referred  for  a  detailed 
derivation.  Note  that  all  of  the  integrations  indicated  are  over  the 
length  of  the  hull  forward  of  the  specified  cross  section. 

The  inertial  forces  and  moments  are  given  by: 


12  ®  /m('n2  +  £ri6  -  zri4)d£ 

(8) 

1 3  =  /m(f)3  -  Cfi5)d? 

(9) 

H  =  J‘(ixn4  -  +  £n6>}d£ 

(10) 

I3  =  -  /m(£  -  x)(n3  -  5n5)d£ 

(11) 

i6  =  /m(5  -  x)(ri2  +  Cn6  -  zn4)d£ 

(12) 

where  ix  =  sectional  moment  of  inertia  about  the  x-axis 
m  “  sectional  mass  per  unit  length 
x  =  x-coordinate  of  the  cross  section 
T  =  z-coordinate  of  the  sectional  centre  of  gravity 
£  =  integration  variable 

r)  j  =  acceleration 
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The  hydrostatic  restoring  forces  and  moments  are  given  by: 


R3  =  -pg/b(n3  -  ^n5)^ 

(13) 

r4  =  gr)^/(pa'55  -  mT)d£ 

(14) 

r5  *  pg/b(?  -  x)(ti3  -  Cn5)d? 

(15) 

with 

R2  =  0  and  Rg  *  0. 

Here 

a  =  submerged  sectional  area 

b  *  sectional  beam  at  waterplane 
g  =  acceleration  of  gravity 
o5  =  sectional  metacentric  height 
rij  =  motion  displacement 
p  =  mass  density  of  water 

The  exciting  forces  and  moments  due  to  the  oncoming  regular  waves 
are  given  by: 

Ej  “  pa{/(fj  +  hj)d£  +  _H  hj(x)}eiwt  j«2,3,4  (16) 

iw 

E5  =  -  pa{/[(^  -  x)(f3  +  h3)  +  ii-  h3Jd^}e^u,t  (17) 

ito 

Eg  =  pa{/{(£  -  x)(f£  +  h2)  +  h2]d^}e^wt  (18) 

iw 

where  U  =  mean  speed  of  ship 

fj  *  Froude-Krilof f  force  (see  Equation  (49)  of  Ref.  2) 
hj  =  diffraction  force  (see  Equation  (50)  of  Ref.  2) 
a  =  wave  amplitude 


by: 


The  hydrodynamic  forces  and  moments  due  to  ship  motions  are  given 


D2  “  “  •r{a22(^2  +  & 6)  +  b22(f,2  +  W 

+  a24*4  +  b24fl4  +  ~7  b22H6  '  Ua22^6}d? 

Uj 

“  lUa22(f,2  +  '  ^7  b22(^2  + 

+  — j(a22^6  +  b22^6)  +  U(a24f)4  "  ^  r  -  x 

a)  a)  * 

D3  “  ~  *^*33^3  "  CK5)  +  b33(in3  " 

"  “7  b33^5  +  Ua33^5^d^  ”  [Ua33(^3  " 

0) 

-  b33(ri3  -  5H5)  -  ^<*33^5  +  b33^5)J  c  =  x 

\  =  "  /{a44^  4+  (t  44+  b  44^  4+  a  $  2+  **  l 

+  b24<n2  +  Cn6)  +  ~  b24n6  -  »«24y 

ur 

-  [Ua24(h2  *  cn6)  -  L.  b24(n2  *  «6) 

*  ^Tu24fi6  *  b2*V  *  U(*44f|4  *  "T  b44f,4)1  r  -  « 
to  W  "» 


(19) 


(20) 


(21) 
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%v/  -V-”  av’a.v 
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,v 


:£•- 


■-  .VjJbwVVW. 


L5  -  /u  -  *M«33<fi3  -  Sn5)  +  b33(n3  - 


J*{Ua33(ri3  -  xri5)  -  b33(ri'3  -  xfij) 

u>2 


-  — ( a33H5  +  b33H5)}d5 
O)2 


D6  =  -/<c  -  x){a22(n2  +  ^fi6)  +  b22^2  +  ^6* 
+  *2u\  +  b24Vd?  '  /{Ua22(fl2  +  XI^6) 


‘  \  b22("2  *  *V  *  ~7<a22^6  *  b22^6^ 


*  »«24n4  -  — r  b24'i4)d? 


where  aj^  =  two-dimensional  sectional  added  mass  coefficient 

bjk  ■  two-dimensional  sectional  damping  coefficient 

t>44*  “  two  dimensional  sectional  viscous  damping  coefficient 

for  roll 


velocity 
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3.  DISCUSSION  OF  RESULTS 


Table  1  summarizes  the  cases  examined  in  the  correlation  study. 

In  the  following  discussion,  only  the  minimum  proportion  of  the  results 
necessary  to  illustrate  the  overall  correlation  characteristics  are 
presented. 

Principal  particulars  for  the  five  model  ships  shown  in  Table  1 
are  given  in  Table  2.,  In  fact,  some  of  the  dimensions  given  are  for  ship 
scale,  like  the  SHIPMO  input,  so  that  SHIPMO  computed  results  would 
correspond  to  published  experimental  ones.  The  reason  for  the  apparent 
inconsistency  in  the  choice  of  the  variables  for  the  abscissas  and 
ordinates  in  Figures  4  through  8  is  to  preserve  the  original  formats  in 
which  the  experimental  results  are  presented  in  the  References,  and  hence 
to  facilitate  the  correlation  of  the  predicted  and  experimental  results. 

Firstly,  the  predicted  moment  and  shearing  force  coefficients  for 
Destroyer  Model  Is*  are  compared  with  the  experimental  results  in 
Figure  4.  Agreement  between  the  predicted  and  measured  vertical  bending 
moments  (VBM)  at  stations  5,  10  and  13  is  good  in  head  seas  (Figure  4a), 
but  poor  in  following  seas  (Figure  4b).  The  predicted  VBM  for  this  model 
tended  to  overestimate  the  experimental  results.  Note  in  Figure  4a  that 
the  magnitudes  of  the  experimental  variations  between  the  pairs  of 
measurements  at  X/L  *  0.6  for  stations  10  and  13,  at  X/L  “  1.2  for 
station  5,  and  at  X/L  =  1.5  for  station  10,  range  from  approximately  30 
to  100%  of  the  lower  values. 

The  predicted  VBM  at  midships  in  oblique  seas  are  compared  with 
experiments  in  Figure  4c.  Correlation  is  somewhat  better  for  quartering 
seas  than  for  bow  seas. 

Figures  4d  and  4e  show  fair  agreement  between  the  predicted  and 
measured  horizontal  bending  moments  (HBM).  Though  the  maximum  HBM  are  not 
clearly  defined  by  the  experiments,  the  theory  tends  to  underestimate  the 
experimental  results.  This  trend  is  opposite  to  that  observed  for  the 
VBM.  Incidently,  it  will  be  noted  in  Figure  4d  that  the  predicted  maximum 
HBM  occurs  in  waves  whose  so-called  effective  wavelength  -  defined  as  the 
wavelength  multiplied  by  the  absolute  value  of  secant  of  the  wave 
direction  -  is  approximately  equal  to  L. 


*  The  numbering  of  ordinate  stations  in  Reference  5  is  slightly  different  from 
that  adopted  in  SHIPMO  in  that  the  station  at  the  forward  perpendicular  (FP) 
is  numbered  1  in  the  former,  while  it  is  to  be  numbered  0  in  the  latter. 

The  numbering  of  ordinate  stations  in  this  memorandum  conforms  to  that  of 
SHIPMO.  Thus,  station  5  is  situated  at  L/4  from  the  FP,  where  L  is  the 
length  between  perpendiculars,  and  station  10  at  midship,  etc.  Incidently, 
Reference  5  contains  the  most  comprehensive  model  test  data  for  wave  loads 
presently  available  in  the  open  literature. 


* 

! 

The  predicted  and  measured  torsional  moments  show  gross 
discrepancies  in  Figures  4f  and  4g.  Further  work  is  clearly  needed  to 
resolve  these.  (Note,  the  scale  division  of  the  ordinate  in  Figures  4f 
and  4g  is  100  times  smaller  than  in  Figures  4a  through  4e.) 

The  comparison  of  the  predicted  vertical  shears  with  experiments 
in  head  seas  is  shown  in  Figure  4h.  The  predicted  and  measured  maxima 
appear  to  be  of  the  same  order  of  magnitude  for  station  5,  where  shear  is 
expected  to  be  greatest.  The  scatter  of  the  experimental  data  is 
^  particularly  evident  in  this  figure. 

Next,  Figure  5  shows  the  comparison  of  the  predicted  and  measured 
VBM  coefficients  at  midships  (station  10)  and  L/10  abaft  midships  (station 
12),  for  Destroyer  Model  II6*  in  head  seas.  (Note  the  experimental 
results  are  presented  by  the  solid  lines  from  Reference  6  rather  than 
points.)  The  predicted  maximum  VBMs  are  seen  to  be  consistently  30  to  40% 

I  greater  than  the  measured  maxima  over  the  Froude  number  (Fn)  range  from 

0  to  0.40.  The  good  agreement  shown  in  Figure  5b  for  station  10  appears 
to  be  fortuitous  judging  from  the  trends  in  the  other  figures  for 
different  values  of  Fn. 

Thirdly,  the  predicted  VBM  and  HBM  coefficients  at  midships  of  a 
|  T-2  tanker  model7'10  are  compared  with  experiments  in  Figure  6. 

•  Overall,  the  predicted  VBM  shows  a  good  correlation  with  experiments, 

though  the  discrepancy  seen  in  Figure  6c  is  appreciable.  On  the  other 
hand,  predicted  HBM  correlates  poorly  with  experiment;  the  theory  tends  to 
underestimate  the  experimental  results  considerably,  and  the  discrepancy 
increases  with  increasing  ship  speed.  A  "weather"  HBM  causes  the  bow  and 
|  stern  to  deflect  toward  the  oncoming  waves;  a  "leeward"  HBM  produces 

‘  deflections  in  the  opposite  direction.  The  analysis,  of  course,  does  not 

differentiate  the  two,  because  the  ship  is  assumed  to  be  executing  the 
simple  harmonic  oscillatory  motions,  linear  and  angular,  about  the 
translating  axes. 

j  The  fourth  case  in  Figure  7  shows  the  comparison  between  the 

predicted  and  measured  VBM  and  HBM  coefficients  at  midships  of  the  Series 
60  (Cg  **  0.80)  hull8'10.  The  agreement  between  the  predicted  and 
measured  VBM  appears  to  be  good  in  nearly  head  seas,  but  poor  in  other 
wave  directions  (Figure  7a).  The  tendency  of  theory  to  underestimate 
experimental  results  in  Figure  7a  is  more  conspicuous  in  Figure  7b,  in 
|  which  the  predicted  HBM  is  shown  to  be  consistently  less  than  the 

experimental  results  in  all  wave  direction. 


*  In  Reference  6,  ordinate  stations  are  numbered  from  the  aft  perpendic¬ 
ular,  AP  (station  0)  to  the  forward  perpendicular,  FP  (station  10)  with 
half-stations  between  integral  stations.  The  station  numbers  cited 
above  correspond  to  the  numbering  system  of  SHIPMO;  so  station  0  is  the 
FP,  and  station  20  the  AP. 
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Finally,  the  predicted  ratios  of  VBM  and  HBM  to  the  wave 
amplitude  at  midships  of  a  cargo  ship,  the  WOLVERINE  STATE9'10,  are 
compared  with  experimental  data  in  Figure  8.  The  agreement  between  the 
predicted  and  measured  VBM  results  is  generally  good  in  head  seas,  but 
becomes  poor  in  other  wave  directions  (Figures  8a  and  8b) •  As  seen  in 
Figures  8c  and  8d,  the  predicted  HBM  results  correlate  poorly  with  the 
experiments. 

Because  wave  loads  on  ships  are  derived  as  a  function  of  the 
ship's  motions  (displacement,  velocity  and  acceleration)  and  the 
hydrodynamic  coefficients  (added  mass,  added  moment  of  inertia,  and 
damping),  the  predicted  wave  loads  will  at  best  be  only  as  reliable  as 
these  variables;  but  when  compounded,  even  relatively  minor  inaccuracies 
in  hydrodynamic  coefficients  and  motions  tend  to  be  amplified  and  result 
in  a  significant  inaccuracy  in  the  predicted  loads.  Comparison  of  the 
trends  of  the  correlations  of  wave  loads  (Figures  4  and  5)  and  of  motions* 
(Figures  A1  and  A2  in  the  appendixt)  pertaining  to  the  two  destroyer 
models  bears  out  this  fact. 

The  results  of  the  present  comparison  show  varied  degrees  of 
agreement  between  theory  and  experiment.  Granted  that,  considering  the 
complexity  of  the  model  tests,  it  will  be  prudent  to  make  a  due  allowance 
for  the  accuracy  of  the  experimental  data  when  assessing  the  reliability 
of  the  theoretical  prediction;  yet,  in  view  of  the  wide  margins  between 
the  predicted  wave  loads  and  the  experimental  results  shown  in  the  present 
comparisons,  the  need  for  further  improvement  of  the  theory  is  evident. 


4.  CONCLUDING  REMARKS 

1.  The  purpose  of  the  present  study  was  to  obtain  information  about 

the  accuracy  of  the  predicted  wave  loads  from  the  extended  SHIPMO  computer 
program,  and  hence  to  determine  its  usefulness  as  a  working  tool  for  ship 
structural  design.  The  wave  loads  for  five  ships  in  regular  waves  were 
calculated  by  SHIPMO  and  compared  with  the  experimental  data.  The  main 
results  were  as  follows: 

(a)  The  agreement  between  the  predicted  and  experimental  results  was 
generally  good  for  the  vertical  bending  moment  in  head  seas,  but 
the  correlation  became  less  satisfactory  for  oblique  and  following 
seas. 


*  Even  though  phase  angles  of  ship  motions  do  not  appear  explicitly  in 
the  equations  for  wave  loads  given  in  Section  2  owing  to  the  fact  that 
motions  are  expressed  in  terms  of  complex  variables,  their  accuracy  is 
equally  important  as  that  of  the  motion  amplitudes  in  the  wave-load 
prediction. 

t  For  a  comprehensive  correlation  study  of  the  predicted  ship  motion  with 
experimental  results,  see  Reference  11. 


11 


(b)  The  predicted  horizontal  bending  moments  were  at  best  in  fair 
agreement  with  experiment. 

(c)  The  predicted  and  experimental  results  showed  glaring  discrep¬ 
ancies  with  respect  to  the  torsional  moment. 

(d)  There  appeared  to  be  no  regularity  in  the  pattern  of  discrepancy 
between  theory  and  experiment. 

No  definite  assessment  could  be  made  about  shears  because  of  lack  of  the 
relevant  experimental  data. 

2.  It  was  encouraging  that  the  predicted  maximum  vertical  moment  in 
head  seas  was  for  the  most  part  in  satisfactory  agreement  with 
experimental  data.  However,  for  the  ultimate  statistical  prediction  of 
the  loads  in  realistic  short-crested  waves,  the  accurate  predictions  of 
the  loads  in  other  headings  are  equally  essential.  Thus,  the  correlation 
between  theory  and  experiment  for  the  loads  in  oblique  waves  would  have  to 
be  further  improved,  before  the  wave  loads  predicted  by  SHIPMO  could  be 
used  with  confidence  for  structural  design  purposes. 

3.  The  present  comparison  illustrates  the  need  for  systematic  model 
tests  with  adequate  descriptions  of  the  uncertainties  in  the  measurements. 
The  published  reports  lacked  this  information.  Consequently,  it  was  not 
possible  to  clearly  define  the  range  of  confidence  of  the  predicted  wave 
loads. 


TABLE  1 :  SUMMARY  OF  CASES  EXAMINED  IN  PRESENT  STUDY 


Range  of  Wave 


Model 

Reference 

Froude  Number 

Heading 

Characteristics 

Destroyer  I 

(5) 

0.21,0.29 

0,  30 
60,  90 
120,150 
180 

0.2  <  X/L 

2.2 

Destroyer  II 

(6) 

0,  0.075 

0.15,0.20 

0.30,0.40 

180 

1.0  <  w/L/g 
<4.0 

T-2  Tanker 

(7,10) 

0  -  0.20 

45,120 

135,150 

180 

X/L  -  1 

Series  60 

CB  *  0.80 

(8,10) 

0.15 

10,  40 
60,  90 
120, 150 
170 

0  <  w/L/g 
<4.5 

WOLVERINE 

STATE 

(9,10) 

0.160,0.214 

0,  30 
60,  90 
120,150 

0  <  X/L 
<1.8 

TABLE  2:  PRINCIPAL  PARTICULARS  OF  MODELS 


Series  60 

Destroyer  I  Destroyer  II  T-2  Tanker  0^0.80 


Wolverine 

State 


«  =  SURGE  tj,  * HEAVE  *  PITCH 

I  3  5 

tjz  =  SWAY  ij4  *  ROLL  ij6  *  YAW 
Figure  1  Definition  of  motions 


V,  «  COMPRESSION  FORCE  V4  «  TORSIONAL  MOMENT 

V,  ■  HORIZONTAL  SHEAR  V5  «  VERTICAL  BENDING 

FORCE  MOMENT 

V.  »  VERTICAL  SHEAR  V6  ■  HORIZONTAL  BENDING 

5  FORCE  MOMENT 


Figure  2  Definition  of  positive  loads 


(e)  Horizontal  Bending  Moment  Coefficient  in  Midships 
(Station  10) 


Figure  4  cont'd 


Figure  4  cont'd 
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(c)  Vertical  Bending  Moment  Coefficient  in 
Head  Seas  Fn  ”  0.30 
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Figure  6  Wave  Loads  at  Midships  of  a  T-2  Tanker  Model 


Effective  Wave  Length  ■  Model  Length 


(c)  Horizontal  Bending  Moment  at  Midships  in  Bow  Seas 


(d)  Horizontal  Bending  Moment  at  Midships  in  Quartering 
Seas 


Figure  8  cont'd 


DESTROYER  MODELS  I  AMD  II 


The  non-dimensionalized  motions  of  destroyer  model  I  predicted  by 
the  SH1PM0  computer  program  are  compared  with  experimental  results  in  Figure 
A1  for  Fn  *  0.21.  The  results  for  Fn  *  0.29  were  similar  in  terms  of  both 
magnitude  and  trend.  In  Figure  Ala,  the  predicted  pitch  shows  a  good 
correlation  with  experiments  in  all  headings,  though  theory  tends  to 
underestimate  the  measured  responses.  The  correlation  of  the  predicted  heave 
motion  with  experimental  results  is  reasonable  for  head  seas,  but  somewhat 
poor  in  oblique  waves;  particularly  in  beam  waves,  the  discrepancy  between 
computation  and  experiment  is  significant  (Figure  Alb).  The  correlation 
between  predicted  roll  and  yaw  with  experimental  data  are  reasonably  good  in 
all  headings  with  the  exception  at  60°  as  seen  in  Figures  Ale  and  d.  In 
general,  the  correlation  between  the  predicted  and  experimental  results  is 
poor  at  long  wavelengths. 


Figure  A2  shows  the  comparison  of  the  predicted  non-dimensionalized 
heave  and  pitch  motions  of  destroyer  model  II  in  head  seas  with  experiment. 
The  results  are  indicated  by  solid  lines  as  done  in  Reference  6.  Except  for 
low  and  high  frequency  tails,  the  correlation  between  the  predicted  and 
experimental  results  is  reasonably  good  over  the  Fn  range  0  to  0.40. 


YAW  AMPLITUDE  ROLL  AMPLITUDE 

WAVE  SLOPE  ~  WAVE  SLPOE 


Roll 


(d)  Yaw 


Figure  A.l  cont'd 
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LEGEND:  7j3  -  HEAVE  AMPLITUDE  /  WAVE  AMPLITUDE 

TJ,  -  PITCH  ANGLE  /  WAVE  SLOPE 

EXPERIMENT  - 

SHIPMO  - 


Fig.  A2.  Motions  of  Destroyer  Model  II  in  Regular  Waves 
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